High temperature superconductivity from realistic Coulomb and Frohlich interactions 
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In the last years ample experimental evidence has shown that charge carriers in high-temperature 
superconductors are strongly correlated but also coupled with lattice vibrations (phonons), signaling 
that the true origin of high-T c superconductivity can only be found in a proper combination of 
Coulomb and electron-phonon interactions. On this basis, we propose and study a model for high-T c 
superconductivity, which accounts for realistic Coulomb repulsion, strong electron-phonon (Frohlich) 
interaction and residual on-site (Hubbard U) correlations without any ad-hoc assumptions on their 
relative strength and interaction range. In the framework of this model, which exhibits a phase 
transition to a superconducting state with a critical temperature T c well in excess of 100K, we 
emphasize the role of U as the driving parameter for a BEC/BCS crossover. Our model lays 
a microscopic foundation for the polaron-bipolaron theory of superconductivity. We argue that 
the high-T c phenomenon originates in competing Coulomb and Frohlich interactions beyond the 
conventional BCS description. 

PACS numbers: 74.20.-z, 74.72.-h, 71.38.-k 



I. INTRODUCTION 

Unconventional symmetries of the order parameter al- 
lowed some researchers to maintain that a purely re- 
pulsive interaction between electrons (Hubbard U) ac- 
counts for superconductivity without phonons in a num- 
ber of high-temperature superconductors . However, re- 
cent analytical 2 and numerical (Monte-Carlo) 3,4 studies 
shed doubts on the possibility of high temperature su- 
perconductivity from repulsive interactions only. 

Also a growing number of experimental and theoreti- 
cal results suggest that strong electron correlations and 
significant electron-phonon interaction (EPI) are the un- 
avoidable features for a microscopic theory of high T c 
superconductivity 5,6. In particular the doping depen- 
dent oxygen-isotope effects on the critical temperature 
T c and on the in-plane supercarrier mass (ref. ), pro- 
vide direct evidence for a significant EPI and bipolaronic 
carriers 1 5 in high-temperature cuprate superconductors. 
Angle-resolved photoemission spectra (ARPES) 14,15 pro- 
vide further evidence for the strong EPI apparently with 
c-axis-polarised optical phonons 1 ''. Some theoretical 
models show that detailed understanding of ARPES re- 
quires EPI 15 ' 17-20 and the lattice disorder 21,22 to be taken 
into account along with strong correlations. These results 
as well as neutron scattering > , tunneling -0-- ' , pump- 
probe 21 ^, earlier 2 " and more recent optical specroscopies 
unambiguously show that lattice vibrations play a signif- 
icant but unconventional role in high-temperature super- 
conductivity. 

Since first proposed , much attention has been paid 
to the strong EPI as a mechanism of superconductiv- 
ity providing effective on-site and inter-site attractions 
between small polarons (electrons dressed by a cloud of 
phonons) 52 . In the framework of negative Hubbard-t/ 



and extended negative Hubbard-U models, the strong 
electron-phonon coupling results in a bound state of two 
polarons that condense with a Bosc-Einstcin critical tem- 
perature strictly related to the mobility of the pairs 53 . 
However, the failure of these models in predicting a high 
critical temperature, due to localization of pairs in the 
strong coupling regime with some particular (Holstein) 
EPI, led to a better understanding of more realistic EPIs 
as the competing interactions with respect to Coulomb 
repulsion. 

Analytical and numerical calculations : 1, 55 clarified 
that EPI with high-frequency optical phonons in ionic 
solids remains poorly screened signaling the presence of 
long-range (Frohlich) EPI at any density of polarons with 
a remarkable reduction of the polaron effective mass. 
Consistently, studies on the so called "Frohlich- Coulomb" 
model (FCM) 36,37 , in which strong long-range EPI and 
long-range Coulomb repulsion are treated on equal foot- 
ing, predict light polarons and bipolarons (bound state 
of two polarons) in cuprates with a remarkably high su- 
perconducting critical temperature in the range in which 
all the interactions are strong compared with the kinetic 
energy of carriers. The interpretation of the optical spec- 
tra of high-T c materials as the polaron absorption 5 ' 29,38 
strengthens the view 36 that the Frohlich EPI is impor- 
tant in those compounds. 

In most analytical and numerical models of high- 
temperature superconductivity, proposed so far, one or 
both Coulomb and electron-phonon interactions have 
been introduced as input parameters not directly re- 
lated to the material. Different from those studies an 
analytical multi-polaron model of high-temperature su- 
perconductivity in highly polarisable ionic lattices has 
been recently proposed 39 and numerically studied (for 
two-particle states) 1 " with generic (bare) Coulomb and 
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Frohlich interactions avoiding any ad-hoc assumptions 
on their range and relative magnitude. It has been 
shown that the generic Hamiltonian comprising any- 
range Coulomb repulsion and the Frohlich EPI can be 
reduced to a short-range t— J p model at very large lattice 
dielectric constant, eo — > oo, for the moderate and strong 
EPI. In this limit the bare static Coulomb repulsion and 
EPI negate each other giving rise to a novel physics de- 
scribed by the polaronic t—J p model with a short-range 
polaronic spin-exchange J p of phononic origin . 

The cancelation of the bare Coulomb repulsion by the 
Frohlich EPI is accurate up to 1/eo corrections. At finite 
eo a residual on-site repulsion of polarons, U, could be 
substantial if the size of the Wannicr (atomic) orbitals 
is small enough. Here we study the effect of this on- 
site repulsion on the ground state of the extended t-J p -U 
model accounting for all essential correlations in high- 
temperature superconductors. It is worth emphasizing 
that the effect of the on-site U does not follows as a mere 
generalization of the t- J p model. The residual Hubbard 
U in fact leads not only to the suppression of on-site 
pairs but also to the reduction of the exchange interaction 
and to the Bose-Einstein condensation (BEC) to BCS 
(BEC/BCS) crossover. 



II. BARE HAMILTONIAN 

Keeping major terms in both interactions, diagonal 
with respect to sites, yield our generic Hamiltonian in 
the site representation, 



H = - ^{TijS^ + ti6ij)4cj + 



h3 



■ / ■ 



huj^hi [tii(q)dq + H.c] + K 



ph • 



(i) 



Here T,,, is the bare hopping integral, if m ^ n, or the site 
energy, if m = n, ^ is the chemical potential, i = (m, a) 
and j = (n, a') include both site (m, n) and spin (a, a') 
quantum numbers, Cj, d q are electron and phonon opera- 
tors respectively, hi = c\ci is a site occupation operator, 
a = e 2 /47re o (e o ~ 8.85 x 10~ 12 F/m is the vacuum 
permittivity), and H p h = X)q^ w q(^q^q + 1/2) with the 
phonon frequency w q . 

The EPI matrix element is 



Ui(q) = (2A^)- 1 / 2 7 (q)exp(iq- m) 



(2) 



with the dimensionlcss EPI coupling, 7(q) (N is the num- 
ber of unit cells ). Deriving the generic Hamiltonian in 
the site representation ! " we approximate the Wannier or- 
bitals as the delta-functions, which is justified as long as 
the characteristic wave-length of doped carriers signif- 
icantly exceeds the orbital size clq. A singular on-site 
(m = n) Coulomb repulsion of two carriers with the 
opposite spins (the Hubbard U) is infinite in this ap- 
proximation. In fact, it should be cut at w a/cooOo as 



indicated by the bar above the sum, J^. Also for mathe- 
matical transparency we consider a single electron band 
dropping the electron band index. 

Quantitative calculations of the EPI matrix elements 
in semiconductors and metals have to be performed nu- 
merically from pseodopotentials. Fortunately one can 
parametrize EPI rather than compute it in many phys- 
ically important cases 41 . EPI in ionic lattices such as 
the cuprates is dominated by coupling with polar optical 
phonons. This dipole interaction is much stronger than 
the deformation potential coupling to acoustic phonons 
and other multipole EPIs. While the EPI matrix ele- 
ments are ill-defined in metals, they are well defined in 
doped insulators, which have their parent dielectric com- 
pounds with well-defined phonon frequencies w q and the 
electron band dispersion. 

To parameterize EPI one can calculate the lowest or- 
der two-particle vertex function comprising the direct 
Coulomb repulsion and a phonon exchange 41 , 



r(q, a,) 



, +| 7 (q)| 2 (^ q ) 2 2%,r!„) . (3) 



Here q = — ki, Cl n = u) n i\ — w n \ are the momentum 
and energy transfer in a scattering process of two carriers 
with the initial momenta and the Matsubara frequencies 
ki^ andwni^, respectively, and 2?(q, fl n ) = — /icj q /[S! 2 + 
(hujq) 2 ] is the propagator of a phonon of frequency w q , 
and Vq is the unit cell volume. In the static limit, 57„ = 
0, Eq.(3) yields the Fourier component of the particle- 
particle interaction as 



T(q,0) 



Ana 



| 7 (q)| 2 So; q . 



(4) 



On the other hand, two static carriers localised on sites 
m and n in the ionic lattice repel each other with the 
Coulomb potential 



(5) 



e |m - n 



where the static dielectric constant, eo accounts for the 
screening by both core electrons and ions. Comparing 
Eq.(4) and Eq.(5) we find 



| 7 (q)| 2 Sw q 



Ana 
KV q 2 



(6) 



at relatively small q < I/a. Here a is the lattice con- 
stant and k = eoe 00 /(eQ — e^) with the high-frequency 
dielectric constant e^. The static dielectric constant eo 
and the high-frequency dielectric constant f m are readily 
measured by putting the parent insulator in a capacitor 
and as the square of the refractive index of the insula- 
tor, respectively. Hence, different from many models of 
high-temperature superconductors proposed so far, our 
generic Hamiltonian with the bare Coulomb and Frohlich 
interactions is defined through the measurable material 
parameters. 
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III. t - J p AND t-J p -U MODELS 



Using the Lang-Firsov (LF) canonical 
transformation 1 -' one can integrate out most of both 
interactions in the transformed Hamiltonian' 50 , 



H = - y fajSaqt +p,Si j )c}c j +H ph + 



i,3 



ijTliTlj , (7) 



since the residual repulsion, v%j, is substantially dimin- 
ished by the large dielectric constant of the polar lattice 
[see Eq.(5)]. Here 



(8) 



is the renormalised hopping integral involving multi- 
phonon transitions with A,; = exp X)q u i(l)^q — H.c. 
and ft, is the chemical potential shifted by the polaron 
binding energy. 

Then using the Schricffcr-Wolf (SW) canonical 
transformation 1 and neglecting the transformed 
Hamiltonian Eq.(7) is reduced to the t — J p Hamiltonian 



as 



39 



n) ( S m ■ S n + \n m n 



(9) 



Here the sum over n ^ m counts each pair once only, 
S m = (1/2)X^ >0 -' cL^^'W is the spin 1/2 opera- 
tor (r arc the Pauli matrices), n m = n m f + w m ^, and 

n m\,i — c int j. Cm T4 are site occupation operators. 

All quantities in the polaronic t-J p Hamiltonian (9) 
are defined through the material parameters, in par- 
ticular the polaron hopping integral, iy = T(m — 
n) cxp[— g 2 (m — n)] with the polaron band-narrowing ex- 
ponent 



3 2 (m) 



2ne 2 



nhwoNVo 



T, 1 - 



cos(q ■ m) 



and 



Jp(m) =T 2 (m)/2 3 2 (m)Sw , 



(10) 



(11) 



It has been proposed that the t-J p Hamiltonian, 
Eq. (9) , has a high-T c superconducting ground state pro- 
tected from clustering' 9 . The polaronic exchange J p is 
attractive for polarons in the singlet channel and repul- 
sive for polarons in the triplet channel. The origin of 
this exchange attraction is illustrated in Fig.l. If two po- 
larons with opposite spins occupy nearest-neighbor sites, 
they can exchange sites without any potential barrier be- 
tween them, which lowers their energy by J p proportional 
to the unrenormalised hopping integral squared. 

Importantly the LF trasfomation Eq.(7) is exact, and 
the SW transformation is accurate for the intermediate 



Figure 1: (Color online) Exchange transfer of two polarons 
with opposite spins between nearest-neighbor sites with no 
potential barrier involved. Horizontal lines illustrate atomic 
levels shifted by the carrier-induced lattice deformation. 



and strong EPI coupling, A > l/y/2z, where A is the 
BCS coupling constant and z is the lattice coordination 
number as discussed in details in Ref. 39. The residual 
repulsion of polarons, in the transformed Hamilto- 
nian, Eq.(7), is small compared with the exchange inter- 
site polaron attraction J p and the short-range bipolaron- 
bipolaron repulsion of about the same magnitude, as long 
as eo 3> a/aJ p . With the typical parameters of the 
cupratcs J p is about 1 eV and a/ a w 4 cV, so that the 
residual inter-site repulsion v mn is small if eo 1 , which 
is well satisfied in all relevant compounds . 

Nevertheless the on-site term in i> mn , Eq.(5), U could 
be substantial, if the size of the Wannier orbitals is small 
enough ao < a. This renormalised U is strongly di- 
minished by the lattice polarization with respect to the 
bare on-site repulsion. We have emphasised in Refs.39,40 
that our model describes carriers doped into the charge- 
transfer Mott-Hubbard (or any polar) insulator, rather 
than the insulator itself, different from the conventional 
Hubbard U or t-J models. The bare Hubbard-?/ on the 
oxygen orbitals (where doped holes reside) in a rigid 
cuprate lattice is of the same order of magnitude as the 
on-site attraction induced by the Frbhlich EPI (w 1 to 2 
eV 44 ), so that the residual Hubbard U could be as large 
as a few hundred meV. We now take it into account in 
the energy of a virtual double occupied state \p) with two 
opposite spins on the same site, 



E n — E„ 



u 



(12) 



Then performing the SW transformation the exchange 
attraction is found as 



Jp(u ' m - n) = ^£ 



k=l 



(2g 2 (m-n))* 
k\(k + u) 



(13) 



where u = U /TujJq. The reduction with respect to 
J p (0,m) is moderate as long as the relative u is less 
than 2g 2 , but becomes substantial for u > 2g 2 , Fig. 2, 
which puts the characteristic bipolaron binding energy 
in the range of a hundred meV comparable with the 
double pseudogap in the cuprates 5 . Importantly J p (u) 
remains large or comparable with the polaron hopping 
integral t = T(a) exp[— g 2 (&)] since the spin exchange of 
the t — J p model , Eq.(ll), does not contain the small 
polaron narrowing exponent cxp(— g 2 ). 
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Figure 2: (Color online) Reduction of the inter-site exchange 
attraction J p (u) / J p (0) by the on-site residual polaron-polaron 
repulsion u = U /Hljq for different values of the polaron mass 
exponent g 2 . 



Hence our extended t-J p (u)-U model including major 
correlations effects reads as follows 

i,j m 

+2 M u > m - n) f S m • S n + -n m n n J (14) 

IV. LOW DENSITY LIMIT AND HIGH T c 

As in Refs. 39,40 we adopt here the strong-coupling ap- 
proach to the multi-polaron problem described by the 
Hamiltonian, Eq.(14), solving first a two-particle prob- 
lem and then projecting the Hamiltonian on the repulsive 
Bose gas of small inter-site bipolarons. Such projection 
allows for a reliable estimate of the superconducting crit- 
ical temperature for low carrier density as long as bipo- 
larons remain small. 

If we neglect the polaronic hopping taking t = 0, then 
the ground and the highest energy states are bipolaronic 
spin-singlet and spin-triplet, respectively, made up of two 
polarons on neighboring sites. The zero-energy states 
[in the nearest-neighbor (NN) approximation] are pairs 
of polarons separated by more than one lattice spacing. 
The on-site bipolaron has energy U > 0. 

For t our exact diagonalization (ED) results on 
finite clusters show that the probability to find NN bipo- 
larons falls as we increase the hopping or the strength 
of the on-site repulsion U as shown in Fig. 3 for a 
100 x 100 square lattice. Consistently, the bipolaron 
size increases but remains on the order of the lattice 
spacing in a wide domain of the parameters (see Fig. 4). 
Importantly, although the small bipolaronic configura- 
tion persists for any values of the hopping at U = 0, 
for U 7^ and large values of t up to a critical value 
t c = UJ p (u)/(2U — 8J p (u)), the presence of a finite on- 
site interaction leads to the crossover from a small to 
a large bipolaronic configuration. Finally, for further in- 
creasing t the system undergoes a phase transition to an 
unbound state at t = t c . The crossover from a small 
to a large bipolaronic configuration is also confirmed by 
the calculation of the bipolaron to polaron effective mass 
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t/J p (u) M p (u) W p (u) 

Figure 3: (Color online) Probability to find two polarons on 
the same site (left panel), on nearest-neighbor sites Pb p (cen- 
tral panel), on more distant sites (right panel) in the ground 
state of the t — J p (u) — U model on a 100 x 100 square lattice 
with different on-site repulsions. 




Small bipolaron regime 
(bound state) 



5 10 _ 15 20 25 30 a 5 1 1 5 2 2 5 30 
U/J p (u) U/J„(U) 

Figure 4: (Color online) Left panel: phase diagram for the 
ground state of the polaronic t-J p (u)-U model on a square 
lattice. Right panel: contourplot of the bipolaron radius r/a 
(a is the lattice constant) for a 100 x 100 square lattice with 
periodic boundary conditions. Different numbers represent 
the value of r/a along the boundaries (dashed lines), empha- 
sizing the increasing of the bipolaron radius as we approach 
the unbound regime. Here r = (|m — n|), m and n being the 
position vectors of the two polarons in the bound state. 



ratio with to** = 2m* in the large bipolaron regime, as 
shown in Fig. 5. 

In the small bipolaron regime, the kinetic energy oper- 
ator in Eq.(14) connects singlet configurations in the first 
and higher orders with respect to the polaronic hopping 
integrals. Taking into account only the singlet bipolaron 
band and discarding all other configurations one can map 
the t — Jp{u) — U Hamiltonian on the hard-core charged 
Bose gas as described in Ref. 39. This gas is superfluid in 
2D and higher dimensions. In particular, its 2D critical 
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W p (u) 

Figure 5: (Color online) Ratio of bipolaron to polaron mass 
in the t — J p (u) — U model on a square lattice. 
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U=0.8eV N. 







0.0 0.2 0.4 0.6 0.8 1.0 



t/J p (u) 

Figure 6: (Color online) The superconducting critical tempe- 
rature of the t — J p (u) — U model on the square lattice for low 
carrier density rib = 0.05/a 2 with frwQ = 0.08eV, g 2 = 1.24 
and J p (0) = leV 39 . 



temperature T c in the dilute limit is given by 40 

rp = 2Trh 2 n b 
c fc s m**lnln(l/n b a 2 ) ' 1 ' 

where rib is the boson density per unit area. 

The occurrence of superconductivity in this regime is 
not controlled by a pairing strength, but by the phase 
coherence among small bipolarons 11 . At low enough den- 
sity the Bose-Einstein condensation (BEC) temperature 
in 3D or its Berezinsky-Kosterlitz-Thouless (BKT) ana- 
log in 2D, Eq.(15) should not significantly depend on the 
bipolaron size as long as it remains small. On the other 
hand increasing U in our model finally results in a bipo- 
laron overlap, where the bipolaron condensation should 
appear in the form of the polaronic Cooper pairs in mo- 
mentum space with a lower critical temperature, rather 
than in real space (BEC-BCS crossover 41 ^' 1 ). Hence, 
we can safely estimate the BEC critical temperature by 
weighting Eq.15 with the probability to find NN polarons 
as ps Pb P (t I J p (u))T c M . As shown in Fig. 6, despite a 
low carrier density, for a physical choice of the parame- 
ters (hw = 0.08eV, g 2 = 1.24 and J p (0) = leV 39 ) the 



critical temperature is found to be well in excess of 100K 
for U = and rapidly decreases with increasing U. 

It is worth noting that, unlike in other theories, the 
strength of the on-site interaction term reduces T£. How- 
ever, we recall that our residual on site interaction U is 
defined as the difference between bare Hubbard U and 
on-site Frohlich EPI therefore at U — 0, when T c is max- 
imized, we have a strong bare on-site interaction with 
U ps 2E„ ~ 1-2 eV. 



V. CONCLUSIONS 

In conclusion, we have introduced and studied the po- 
laronic t-J p (u)-U model, defined through the bare mate- 
rial parameters. The model, being an essential general- 
ization of the t-Jp model 4 ', includes all electron-electron 
and clcctron-phonon correlations providing a microscopic 
explanation of the high-T c phenomenon without any ad- 
hoc approximations. We show that the inclusion of 
the residual on-site interaction U (neglected in the t-J p 
model 39 ' 40 ), drives the system to a BEC/BCS crossover 
that reconciles the polaron-bipolaron theory of supercon- 
ductivity with the observation of a large Fermi surface 
in overdoped cuprate superconductors. We offer an ex- 
planation, on microscopic grounds, of the high-T c phe- 
nomenon as a consequence of competing Coulomb and 
Frohlich interactions in highly polarizable ionic lattices 
beyond the conventional BCS description. 
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